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Abstract 

Q\ \ We have computed the 77 -pseudoscalar octet mass splitting using staggered 

fermions on both dynamical and quenched gauge configurations. We have 
used Wuppertal smeared operators to reduce excited state contributions. We 
compare our results with the theoretical forms predicted by partially quenched 
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chiral perturbation theory in the lowest order. Using lattice volumes of size 
16 3 x 32 with a -1 = 2GeV we obtain results consistent with the physical rf 
mass. We also demonstrate that the flavor singlet piece of the rf mass comes 
from zero modes of the Dirac operator. 

I. INTRODUCTION 



By now there have been dozens of calculations of the masses of most of the light hadrons 
using lattice QCD. However, until recently the rf meson received only scant attention, in 
large part because of the relative difficulty of the calculation: the disconnected contraction 
Ph! which gives the vf propagator its special character is an order of magnitude more expensive to 
compute. A previous study [|l[] made a first study of the rf two-point function using quenched 
Wilson fermions, finding a vertex of the right size to explain the rf mass. Here we improve 
the situation in several ways, by using staggered fermions which have better chiral properties, 
by going closer to the chiral limit, and most importantly by using gauge configurations with 
dynamical fermions. The authors of the earlier study also took a step toward confirming 
the conventional wisdom that the rf receives its special mass from instantons by sorting 
their gauge fields into bins of topological charge. Here we take the further step of examining 
the contribution of topological zero modes of the Dirac operator to the rf propagator. The 
results again confirm the lore that this is the mechanism by which the rf gets its mass. 

The rf meson is one of the more intriguing strongly interacting particles. Although it is 
the lightest flavor singlet pseudoscalar, it is too heavy to be the Goldstone boson of a U(l) 
axial symmetry, as was pointed out by Weinberg long ago. With the emergence of QCD 
it was understood that this U(l) symmetry is actually anomalous, so that one should not 
expect an associated Goldstone boson. A qualitative understanding of the mechanism behind 
the rf mass was provided in the framework of the large N c expansion. When N c — > 00 , 
£7(1) axial symmetry is restored and gives rise to rf — 7r degeneracy. This degeneracy is lifted 
by the presence of virtual quark loops (Fig |l|) in the rf propagator which are suppressed by 
one power of 1/N C . Infinite iteration of these quark antiquark annihilation diagrams gives 
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rise to a geometric series for the Euclidean rf propagator. Defining the correlator (i]' (t)r]' (0)) 
one can write its Fourier transform as 



1 mi 1 2 1 2 1 
W ~ 2 5— 2 H 2~ 2 m o^- 2 m o^- 2 + I 1 ) 

+ TOg £T + TOg p z + Wig £T + nig + 77lg 

where m§ is the average of the square of the octet masses and m\ is the strength of the flavor 
singlet interaction. Summing the geometric series shifts the pole from m\ to m\ + m\. Thus 

faW(o)> 



jr + m§ + mg 

In an £77(3) symmetric world, one would write 

m\, =m 2 + ml (2) 

wig vanishes in the chiral limit while m\ does not. Taking m\ = (Am 2 K + 3m 2 + m 2 )/8 
and substituting the experimentally measured values for m v >, m K , m v and m v , one obtains 
m = 860MeV. 

The large N c approximation also gives the qualitative understanding that this mo is 
linked to instantons. To leading order in 1/N C one finds the Witten-Veneziano formula 

ml = ^ (3, 

J n 

where Nf is the number of flavors, x is f ne topological susceptibility and f n is the pion decay 
constant. In the real world, taking the number of active light flavors to be three, one needs x 
of order (185MeV) 4 to explain the rj'. Lattice QCD calculations of the quenched topological 
susceptibility (e.g. 0) do indeed give results in this range. Thus modulo concerns over 1/N C 
corrections, definitions of topological charge, setting the quenched scale etc., one could claim 
that lattice calculations give a quantitatively correct indirect determination of mo- 



II. DIRECT LATTICE CALCULATION OF M 

Lattice QCD offers the challenge of calculating m\ directly from the rj' propagator with- 
out resort to any argument based on large N c . The quantity of interest is the two point 
function (r] r (x)r)'(y), where the simplest choice for the interpolating field T)'(x) is Q(x)75 Q(x). 
As illustrated in figure 0, there are two types of contractions to take into account: (i) the sin- 
gle quark loop connected diagram which contributes to flavor- singlet and flavor non-singlet 
mesons alike, and (ii) the disconnected diagram which appears only for the flavor singlet rj'. 
As in the previous lattice studies |],|6|-f!, we find it convenient to define the ratio R(t) of 
the disconnected two loop amplitude to the connected one loop amplitude 

Noting the asymptotic behavior 
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(rf(0)rf(t)) Z'exp(-rrtyt) 
(7r(0)7r(t)) Z exp(— m n i) 



(5) 



one sees that in full unquenched QCD the ratio R(t) asymptotes to 1 — B exp(— Amt), where 
B is a constant and Am = m r /—m%. This statement assumes that there are equal numbers of 
dynamical and valence flavors and md yn = m va i. If one has differing numbers of valence and 
dynamical fermions, as we do when using staggered fermions, then one needs to rescale the 
connected diagram of fig § by Nf/N v , and the disconnected diagram by (Nf/N v ) 2 . Therefore 
R(t) takes the form 

R{t) = jf [1 - 5exp(-Amt)] (6) 

On quenched configurations, the absence of closed quark loops means that the basic 
vertex is not allowed to iterate, and equation [I] is truncated after the first two terms. Thus 
in the quenched approximation, there is a double pole in the rj' propagator. This means for 
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the zero spatial momentum state, R(t) rises linearly with the slope ^^t. 
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R(t) = const. + — 2_ t (7) 
2m 8 

We have used staggered fermions, both dynamical and quenched configurations and local 
and Wuppertal smeared operators for the extraction of m v i. We also study the dynamical 
flavor dependence of m v i and have derived the expected theoretical forms for the ratio 
using partially quenched chiral perturbation theory (PQxPT) so as to enable comparison 
of our results with theory. In section 3, we review the basic concepts of PQxPT [§] and 
derive expressions for R(t). In section 4, we describe details pertaining to the parameters of 
the simulation and the Wuppertal smearing procedure. In section 5 we discuss the results 
obtained from our simulation and compare with the theoretical predictions of section 2. 



III. RATIO FROM PQxPT 

Bernard and Golterman [|T0|,|9f have developed a technique for constructing an effective 
chiral theory for quenched and partially quenched QCD. The basic idea in this approach is 
based on the observation that if a scalar quark (q) is added to the QCD Lagrangian, then the 
scalar determinant can cancel the quark determinant. For concreteness, if we assume there 
are N v quarks of masses m^i = 1,2..N V ), then adding k pseudoquarks of masses mj(j = 
l,2..k) such that mj = mi would mean that first k quarks are quenched and the remaining 
N v — k unquenched. The low energy effective chiral theory then describes interactions 
among all kinds of boundstates, including the ordinary pions as well as unphysical states 
containing scalar quarks. The reader is referred to pH|J9|] for discussion regarding the form of 
the Lagrangian and its symmetry properties. For the purposes of this paper, it is sufficient to 
write down the Green's function in momentum space in the basis of the states corresponding 
to uu, dd and their pseudoquark counterparts: 

_ dij€j m\ 

ij p 2 + M? (p 2 + M?)(p 2 + Mj)F(p 2 ) 1 J 
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where 



d=k+lP + M d 



(9) 



M 2 is the square of the meson mass composed of quark flavor i and 

+1 1 < i < N v 
-1 N v + l<i< N v + k 



(10) 



In the above, the term proportional to a, the momentum dependent self-interaction of 
r]' has been neglected. It can be reinstated any time by the substitution Wq — > m\ + a 2 p 2 . 
When rrii = rnj the first term is simply the neutral meson propagator in the absence of 
flavor singlet interactions. The second term is obtained from summing the infinite geometric 
series obtained by iterating the flavor singlet interaction. When simulating with staggered 
fermions on Nf = 2 dynamical configurations, the variables N v and k take the values 4 and 
2 respectively. Specifically, when all the valence flavors have the same mass m va i and all the 
dynamical flavors have mass md yn but m va i ^ rridyn, (in a typical simulation this situation 
would arise if one were simulating with staggered fermions on dynamical configurations) the 
Green's function above takes the form, 
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ml (p 2 + M 2 ) 



P 



M? (p 2 + M 2 ) 2 {p 2 + M 2 + N f m 2 ) 



(11) 



It is straightforward to go over to configuration space, project each term on to zero spatial 
momentum state and then take the ratio of the second term to the first termQ. Doing so we 
obtain 



R(t) 



Nv 
Nf 



[At - Bexp(-Amt) + C] 



where 
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{m 2 d -ml){m 2 d - m 2 ,) 
(ml - m 2 ,)2m 8 



(12) 



(13) 



B 



m 8 (m^ - m 2 ,) 2 
m ri i(ml — m 2 ,) 2 



(14) 
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[ml 



m 2 ,) 



2ml 
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(ml - m 2 ,) 



+ 1 



(15) 



Mf = M: 



3 ' 



m„ 



M] and 



(16) 



All the formulae have been derived in the infinite volume limit 
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ni 



m 2 d + N f m 2 . (17) 



As a check we note that eqn [12] contains the expression in eqn ^| and eqn as appropriate 
limits. When m 8 = one obtains eqn |6] with B = m s /m v r, and in the limit — > oo, 
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Nf — > one obtains eqn [7| with the intercept equal to We have ratio data on quenched 
and Nf = 2 configurations. Our ratio data on Nf = 2 configurations can be further divided 
into m va i = rridyn and m va i ^ md yn - Quenched and m va i = rrid yn data can be subjected to 
both one and two parameter fits thus enabling comparison between the predictions of the 
theory and "experiment" . Likewise, a similar comparison is possible from both one and four 
parameter fit to m va i ^ m dyn data. 



IV. SIMULATION DETAILS 



A. Ensemble 



The propagators required for computing the disconnected and the connected amplitude 
were computed on configurations of size 16 3 x 32. The statistical ensemble used and the 
valence quark masses at which the simulation was performed are shown in table |. The 
dynamical configurations were borrowed from Columbia while the quenched configurations 
were generated on the Cray T3D at Ohio Supercomputer Center (OSC). By design, the 
ensembles listed in table | have the same inverse lattice spacing of about 2GeV obtained 
from m p . The values of m va i on the quenched configurations were chosen 10% higher than 
those on the dynamical configurations, in order to bring the corresponding values of m 8 
into more precise agreement. The staggered propagators were computed by invoking the 
conjugate gradient method using sets of processors ranging from 16 to 64 nodes of the 
Cray T3D at OSC and Los Alamos's Advanced Computer Laboratory with a sustained 
performance of 45 Mflops per node. 



B. Propagators 



The pseudoscalar operator that creates or destroys an r/ in the staggered formalism is 
Q"i§®IQ- It is a distance 4 operator in which the quark and the antiquark sit at the opposite 
corners of the hypercube. We put in explicit links connecting the quark and the antiquark 
across the edges of the hypercube and averaged over all the 24 possible paths to make the 
operator gauge invariant. As was mentioned earlier, the two point correlation function of 
this operator yields both the connected one loop and the disconnected two loop amplitude 
of fig ||. Two propagators need to be computed for calculating the one loop amplitude. 
One propagator was calculated using a delta function source with a local random phase 
7] x = exp(i6(x)) at each spatial site and for each color on time slices t = 1 and t = 2. When 
averaged over all noisy samples, the noise rj x satisfies (rj x vl) = *W Translating the noisy 
source at each site (made gauge invariant by putting in links) a hypercubic distance 4 and 
putting in the phase appropriate for the operator, the second propagator was calculated. 
Transporting the source and calculating the anti-propagator separately is the price paid due 
to the non locality of the staggered flavor singlet operator. For each configuration we took 
two noise samples with the lattice doubled in the time direction. 
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Very few lattice calculations of exist because of the difficulty involved in getting good 
signals for the disconnected amplitude using reasonable computer time. We have addressed 
this problem by using a noisy source like the one used for the connected amplitude but placed 
on all sites of the lattice. Then we solve (p + m)0 = 77 and estimate the quark propagator 
as G xy = (m(f) x <jy y ) with \x — y\ even. The distance 4 staggered flavor singlet operator that 
we use satisfies this criterion. At this point it should be mentioned that there also exists a 
distance 3 staggered flavor singlet operator for which the the appropriate estimator of G xy 
is (r)l<f)y). This channel did not yield good signals for the the pseudofermion propagators, 
and we did chose not to examine it further. In any case, (m<j) x <j) y ) is a better estimator than 
(Vxfiy) srnce the fluctuations in the former go as (xx) 2 while those for the latter are of the 
order (xx)/ m vai where (xx) is the condensate. For the range of light quark masses chosen 
in this simulation, (xx) ues between and 1, thus making (xx) 2 significantly less than 

m-val 

We used 16 noise samples per configuration, again on a doubled lattice. 



C. Wuppertal Smearing 

Fig [| shows the ratio R(t) on quenched and dynamical configuration (Nf = 2, m dyn = 
0.01) using local interpolating fields. Since the disconnected data is noisy, the ratio data 
is only useful for the first few time slices where the contribution due to excited states is 
non- negligible. Smeared operators reduce the contributions from excited states and enable 
a more reliable extraction of Am from the available data. 

We used the Wuppertal smearing technique Jl lj in which one obtains an exponentially 
decaying bound state wavefunction from solving 

(-V 2 [U} + K 2 )<f)(x) = 6 x , (18) 

where k 2 is a parameter that can be tuned to control the spread of the wavefunction. On 
the lattice, for staggered fermions and to make the smearing procedure gauge invariant, the 
operator V 2 takes the form 

1 3 

V 2 [U] -> - E U»(x) U„(x + n) 8 X , >X+ ^ + U£(x - fj) Ul(x - 2/i)4v- 2/ x - 28 x , fX (19) 

We experimented with different values of k 2 and determined that the critical value oc- 
curred near k 2 = —0.64. The values of k, 2 that are relevant to this study are those near 
the critical value since this corresponds to a maximum spread in the wavefunction with- 
out losing the exponential behavior. The form of the wavefunction obtained on one of the 
Nf = 2(m ( iyn = 0.01) configurations is shown in fig |j. 

For this study, one propagator of the connected correlator was calculated with the source 
smeared to a fixed radius corresponding to k, 2 = —0.6 and a point-like sink. This was 
tied to five different anti-propagators with point-like source but smeared at the sink end 
corresponding to different values of k 2 (see fig |5|). Thus the smeared valence propagators 
that we have calculated correspond to the set of correlators shown in table [TI]. 

A typical effective mass plot obtained on Nf = 2, md yn = 0.01 configuration at m va i = 
0.01 is shown in fig ^|for the correlator LL. Since decrease in k 2 increases the spread in the 
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wavefunction, among the 5 correlators shown in table [Tj], C5 is expected to perform the best 
which is evident by comparing fig ||] and fig [7[ We also defined the correlator 

(cos9Q(x)Q(x) + sm9Q w4 {x)Q{x)) (20) 

which amounts to taking the linear combination 

(cos 2 9LL + sm 2 9C 5 + 2cos9 sin0Ci) (21) 

We determined the angle 9 for which the value of the effective mass on t = 1 was the lowest 
and this corresponded to 9 = .7977T. We denote this linear combination of correlators as 
LLC5 and it is clear from fig |7| that it is slightly better than C5. 

To calculate the corresponding disconnected contributions of the correlators listed in 
U the sink end of the pseudofermion propagator (</>') was smeared to 4 different radii 
corresponding to k 2 = —0.5, —0.53, —0.56 and —0.6 (see fig ||). Accordingly, the quark 
propagator was obtained from the estimator G' xy = (m(j) x 4> y ). 

Fig [| compares the ratio plot with and without smearing on dynamical configurations. 
In the initial time slices both the data are different but they begin to coincide after a few 
time slices as they should. 



V. RESULTS 

In the plot of figure |9|, the different curves represent the fit to the ratio data obtained 
from the various configurations listed in table |. For Nf = 2 the points shown correspond 
to m va i = rridyn while for the quenched data the points correspond to m va \ = 0.011. The 
quenched data is fit linearly according to equation [7| with both the intercept and the slope 
as free parameters. The dynamical data is fit according to equation ^| with both and 
B as free parameters. A remarkable observation that is to be gleaned from this graph is 
that the ratio data is distinctly different for different numbers of dynamical flavors and in 
accordance with the predicted theoretical form. If nothing else, this observable is clearly 
able to distinguish quenched from dynamical configurations. 

Our values of m 8 calculated on quenched and Nf = 2 dynamical configurations are 



plotted as a function of the valence quark mass in fig [L0[ We note that at this finite lattice 
spacing, m§ does not quite vanish in the chiral limit. This is a consequence of our taking the 
flavor singlet (75 ® I) pion as opposed to the special (flavor T 5 ) staggered pion. Of course 
even for our pion, the intercept should vanish in the continuum limit. 

From the linear fit to the quenched data for m va i = 0.011, 0.022 and 0.033, we extract 
itLq from the slope. A plot of m\ against m va i and the extrapolation to the zero quark mass 
is shown in fig [11]. m\ increases with decreasing quark mass and does not vanish in the chiral 
limit. Smearing has lowered the value of the intercept obtained in the zero quark mass limit, 
bringing it closer to the real world value. 

Two parameter fits, according to eqn of the dynamical m va i = m dyn data yields Am 
and hence m rj > . Extrapolating the values of m v > so extracted, we find that it does not vanish 
in the chiral limit (fig |l^). As seen in quenched case, the values of m v > both at finite and 
zero quark mass due to the correlator LLC5 are lower than the corresponding points due to 
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the local correlator LL. Errors quoted on the ratio data and all the mass values are obtained 
by doing a single elimination jackknife. 

We could not obtain stable 4 parameter fits to our m va i ^ m dyn data. The known theo- 
retical forms for the ratio obtained from lowest order PQxPT were used for one parameter 
fit of the m va i 7^ mdyn an d m va i = md yn data. Fig [13| shows the theoretical fit obtained for 
m va i = 0.02 on Nf = 2 , mdyn = 0.01 configuration using smeared operators. One param- 
eter fit to the data obtained with local operators gave reasonable \ 2 only if the fit range 
began from t > 3. This suggests that the local data is heavily contaminated with excited 
states since the formulae derived in section 2 hold for asymptotic times. Doing a similar job 
for m va i = 0.03 data and extrapolating the value of m v > obtained from both to the chiral 
limit we obtain m v >(Nf = 3) = 876 ± 16MeV which is remarkably consistent with the value 
obtained from fully quenched and dynamical data (table [m]) . 

A one parameter fit to m va i < m^ n data did not yield reasonable x 2 - However, for the 
m vai = m d yn data, it was possible to extract the deviation from the lowest order PQxPT 
as determined by the observable Z'/Z, the ratio of the residues for creating 77' and 77. 



Z' 
~Z 



1 <o| 


i/(0) 


1) 


2 


1(0 


1*7(0) 


v) I 


2 



(22) 



In PQxPt, this ratio is 1 based on the assumption f v > = f v . Therefore, in PQxPT as shown 
in section 2, the parameter B in eqn ^| equals instead of ^ z - From the 2 parameter 

fit to the m va i = m^yn data one can then determine the ratio Z'/Z. The values we extract 



for Z'/Z are plotted as a function of m va i in fig [L4| both for smeared and local data. It can 
be seen that the data prefer to be 20-30% above unity which is typically the size of higher 
order chiral and 0(1/N C ) corrections. 

In the analysis so far, we have neglected the effect of momentum dependence of m$, ie. 
we had been working in a theory with a, the coefficient of the kinetic energy term of rj' in 
the chiral Lagrangian ||, set to zero, a is expected to be a small quantity, contributing only 



at the next to leading order in a combined expansion in 1/N C and quark mass [10j . Including 
a simply shifts the strength of the flavor singlet interaction vertex from m\ to m\ + a 2 m\. 
While the form of R(t) remains the same as before for all the three cases, the coefficients of 
the time dependent terms and the constants become functions of two unknown parameters, 
ml and a. For the quenched case, we obtain 



(mi — ami) t (mi + ami) 

R(t) = ^— ° ZJ- + - — 23 

2m 8 2 ml 



It is clear from the above formula that the value of m\ is shifted by a small amount at finite 
quark mass and indeed this is borne out true by the quenched ratio data when fit to the 
above form. We see 10%-15% downward shift in the the values of m\ obtained from the 
slope and the intercept of the fit at non zero quark mass. When extrapolated to the chiral 
limit, the value of ??1q is left unaffected (within quoted errors in table |T|) by the presence 
of the momentum dependent interaction term. The statistical errors associated with the 
extracted values of a are rather large. Our results for a for each quark mass is shown in 
table fV]. 

On dynamical configurations, the presence of a non zero a is reflected in the expression 
for m v >. It now takes the form 
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The natural quantity that is extracted from fitting the dynamical m va i = md yn data with 
equation ^| is m v ' which cannot be used to determine both and a independently. 



VI. FERMIONIC ZEROMODES 



As reassuring as it may be that our lattice results correctly reproduce the physical rj' 
mass, the calculation itself does not illuminate the precise mechanism by which the rf gets 
its special mass. The conventional lore, as quantified by the Witten-Veneziano formula, 
associates m with the fluctuations in the topological charge. In a previous study, Kura- 
mashi et al. tried to verify this relation by calculating the topological susceptiblity of their 
gauge configurations and using it to calculate m . They obtain a value that is higher than 
the experimantal number. Here we chose to illuminate the question by focusing directly 
on the fermionic zero-modes which are associated with the topological charge of a gauge 
configuration via the index theorem. 

In the continuum, integrating the anomaly equation gives the relation 



where Sf is the fermion propagator. Resolving the propagator in a sum over eigenmodes of 
p, and noting that all modes with nonzero eigenvalue come in conjugate pairs, one gets the 
index theorem: 



Here n + {nJ) is the number of positive (negative) chirality zeromodes. If the 7/ is particularly 
sensitive to the topological charge of a configuration, then one must also expect it to be 
sensitive to the presence or absence of fermionic zeromodes. 

On the lattice, both sides of the index equation are slightly distorted. The flavor singlet 
75 does not exactly anticommute with p, so the trace of 75SV cannot be collapsed into a 
trace in the zeromode sector alone. Further, there are no exact zeromodes, nor is there an 
exact definition of the topological charge. Nevertheless, if we are sufficiently close to the 
continuum, something like the conventional picture should obtain. 

We have taken a look [|T3J by constructing the lowest few eigenmodes of p on both the 
quenched and dynamical (Nf = 2, m dyn = .01) ensembles of gauge fields. We use the method 
of subspace iterations as described in ref. |T2[]. This algorithm obtains the eigenmodes of p 2 
by minimization of the Ritz functional 



Here X e is the rectangular matrix of m eigenvectors. The subscript indicates we only 
need the solution on the even sites; full eigenmodes of p can then be easily reconstructed 
on the odd sites. We found it necessary to modify the algorithm slightly, periodically 
reorthogonalizing the columns of the block vector X e . With the modes themselves in hand 



Q = mti (75 S F ) 



(25) 



Q = n + — n_. 



(26) 




(27) 



one can then construct approximate propagators and ask how well these small eigenvalue 
modes reproduce the full answers we obtain by conventional means. 

In figure [15] we show the contribution to the topological charge coming from the eigen- 



modes in one typical dynamical configuration. In terms of the eigenmodes, equation |25| takes 
the form 

Q = m£4^- (28) 

As can be seen, after the first few modes, the sum quickly saturates to the full answer 
obtained by averaging over many copies of pseudofermion noise. 

That this behavior is typical of the whole ensemble can be seen in the scatter plot of 



figure which plots the eigenvalue, (75) pairs collected on all configurations. The largest 
contributions to the trace evidently come from the smallest eigenmodes. We conclude that 
even at this finite lattice spacing, the index theorem is perfectly recognizable. 

It is only natural to extend the analysis further to calculate the 7/ disconnected corre- 
lator, which is proportional to the fluctuations in the topological charge and is given by 
(tr(7 5 G(x, x))tr(7 5 G(y, y)))- The result of calculating this using the available eigenmodes 
is shown in figure It is evident that the first few modes give essentially the full answer. 

One may wonder if the chiral modes may be used to successfully reproduce other hadronic 
correlators as well. In terms of the eigenmodes, the approximate propagator is of the form 

E «£M(£) (29) 



Figure [18| answers this question for the case of the Goldstone pion. Even with the full 32 
modes, the pion propagator is off by a large factor both in its overall normalization, and in 
its mass parameter. Evidently the 77' is special in that it is exquisitely sensitive to the lowest 
handful of modes. 



VII. CONCLUSIONS 

We have computed both the quenched flavor singlet vertex m and the full dynamical 
7/ mass, finding agreement with the experimental numbers. Since the statistical errors 
are relatively large we have not attempted to pin down the systematic errors, e.g. from 
finite lattice spacing and from the neglect of rj-rj' mixing. Our results both confirm our 
conventional understanding of the rf in QCD and show that lattice QCD can be a useful 
tool in this difficult sector. We have also shown that unlike other hadron correlators, the 
7/ propagator is particularly sensitive to the presence of fermionic zeromodes. This result 
supports the connection of the rj' to topology. 
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TABLES 



TABLE I. The Statistical Ensemble 



A T 

AT) 


m dyn 




^ samp 


«W 





OO 


6.0 


83 


0.011 










0.022 










0.033 


2 


0.01 


5.7 


79 


0.01 










0.02 










0.03 


2 


0.015 


5.7 


50 


0.01 










0.015 


2 


0.025 


5.7 


34 


0.01 










0.025 



TABLE II. Smeared valence operators. The spin/flavor is 75 <g) I in each case. 





Contraction 


k 2 at source 


K 2 at sink 


LL 


Q(x)Q(x) 


Q(0)Q(0) 


no smearing 


no smearing 


Ci 


g(x)Q(x) 


Q w4 (o)Q(o) 


-0.6 


no smearing 


c 2 


Q w i( x )Q( x ) 


Q w4 (0)Q(0) 


-0.6 


-0.5 


c 3 


Qw2(x)Q(x) 


Q w4 (0)Q(0) 


-0.6 


-0.53 


c 4 


Qw3( x )Q( x ) 


Q w4 (0)Q(0) 


-0.6 


-0.56 


c 5 


Q w4 (x)Q(x) 


Q w4 (0)Q(0) 


-0.6 


-0.6 



TABLE III. m 7) >{N f = 3) in the chiral limit. 




Quenched 


Dynamical (m dyn = m val ) 




(MeV) 


(MeV) 


LL 


1156(95) 


974(133) 


LLC5 


891(101) 


780(187) 


TABLE IV. a from the quenched ratio data. 


m va i a 


LL 


LLC 5 


0.011 


-1.7496(2321) 


-0.2988(2091) 


0.022 


-0.8475(1301) 


-0.1431(1195) 


0.033 


-0.4645(910) 


-0.0399(702) 
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FIG. 1. Virtual quark loops in the rj' propagator 
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FIG. 2. Contributions to 7/ correlator 
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FIG. 4. Wavefunction from Wuppertal smearing. 
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FIG. 6. Effective mass plot for mg with LL correlator. 
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FIG. 7. Effective mass plot with C5 and LLC5 correlators. 
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FIG. 8. Dynamical ratio from local and smeared operators. 
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FIG. 9. Flavor dependence of R(t). 
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FIG. 14. Z'jZ vs. m val . 
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FIG. 15. Topological charge on a typical dynamical configuration 
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FIG. 16. Above: ^5) on the dynamical ensemble. Below: Fluctuations in 
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FIG. 17. Below: Comparison of the two loop amplitudes. Above: Ratio of the full two-loop 
amplitude to that calculated with modes. 
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FIG. 18. Pion propagator on dynamical ensemble at ma = 0.01. 
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